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Abstract 


We examine gravitational wave memory in the case where sources and detector are in an ex¬ 
panding cosmology. For simplicity, we treat the case where the cosmology is de Sitter spacetime, 
and discuss the possibility of generalizing our results to the case of a more realistic cosmology. We 
find results very similar to those of gravitational wave memory in an asymptotically flat spacetime, 
but with the magnitude of the effect multiplied by a redshift factor. 
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INTRODUCTION 


Gravitational wave memory, a permanent displacement of the gravitational wave detector 


after the wave has passed, has been known since the work of Zel’dovich and Polnarev 


extended to the fnll non 


inear theory of general relativity by Christodoulon , and treated 


Q. 


by several anthors [sHlSj. As is nsnal in the treatment of isolated systems, all these works 


consider asymptotically flat spacetimes. However, we do not live in an asymptotically flat 
spacetime, bnt rather in an expanding nniverse. For sonrces of gravitational waves whose 
distance from the detector is small compared to the Hnbble radius, modeling the system as 
an asymptotically flat spacetime should be sufficient. However, some of the most powerful 
sources of gravitational waves (e.g. the collision of two supermassive black holes following 
the merger of their two host galaxies) are at cosmological distances where the asymptotically 
flat treatment is not sufficient. 


In this paper we will treat gravitational wave memory in an expanding universe. To 
avoid the complications of the full nonlinear Einstein equations, our treatment will use 
perturbation theory. There is a well developed theory of cosmological perturbations (see 
e.g. the textbook treatment in j^). However, this standard cosmological perturbation 
theory uses metric perturbations, and we have found llj] that the properties of gravitational 
memory are made more clear when using a manifestly gauge invariant perturbation theory 
based on the Weyl tensor. Cosmological perturbation theory using the Weyl tensor was 


developed by Hawking [17|. We will use a treatment similar to that of [l7j, but also. 


using the conformal flatness of Friedman-Lemaitre-Robertson-Walker (FLRW) spacetimes. 


a treatment that draws heavily on the techniques used in 


ll|. 


Cosmological perturbations depend on the equation of state of the matter. The universe, 
both at the current time and at any previous times from which a realistic source of gravi¬ 
tational wave memory could come, is dominated by dust and a cosmological constant. For 
simplicity, in this treatment we will treat only the case of a cosmological constant, leav¬ 
ing the more general case of dust and a cosmological constant for subsequent work. Thus 
this work treats gravitational waves in an expanding de Sitter spacetime. The perturbation 
equations are developed in section H, the cosmological memory effect is calculated in section 
HI, and the implications of the results are discussed in section IV. 
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EQUATIONS OF MOTION 


From the Bianchi identity 'W\eRap \'^5 = 0 we have 

( 1 ) 

where Saji = Raji — ^R 9 a /3 and R = g°‘^Ra^. Using the Einstein held equation with cosmo¬ 
logical constant 

Raji 2R9otl3 R R^afi ‘^^R'aji ( 2 ) 

we hnd that eqn. o becomes 

eCaji'yS = 87rV[^X5]y3 (3) 

Where - \Tgaii and T = g'^^Tap. 

Both the Weyl tensor, and Tap vanish in de Sitter spacetime. It then follows that when we 
perturb eqn. ([3]) from a de Sitter background, the perturbed equation takes the same form 
with the Weyl tensor and stress-energy replaced by their (gauge invariant) perturbations 
and the metric and derivative operator replaced with their background values. We will 
rewrite this perturbed equation in a convenient form making use of the conformal hatness 
of de Sitter spacetime. Recall that the line element in a spatially hat Friedman-Lemaitre- 
Robertson-Walker (FLRW) spacetime takes the form 

ds^ =—dt^ + a^{t){dx^ + dy‘^ + dz^) . (4) 


Then introducing the usual conformal time g hj g = J dt/a, we hnd that the line element 
takes the form 


ds^ = 


—dg^ + dx^ -f dy^ + dz^ 


(5) 


That is, the de Sitter metric takes the form gap = (Rgap where gap is the Minkowski metric 
with Cartesian coordinates {g,x,y,z). It then follows that the perturbed eqn. ([3]) takes the 
form 


(a-^a 


aPyS 


= Svr 


ad[^Xs]p + Xpi^ds]a + gp[jXs]xd^a 


( 6 ) 


Here da is the coordinate derivative operator with respect to the Cartesian coordinates 
{g, X, y, z). Also here and in what follows we use the convention that indicies are raised and 
lowered with the Minkowski metric gap- 
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Following the method of Hj] we now decompose all quantities in terms of spatial tensors 
as follows, using latin letters for spatial indicies: 



Eab — ^ Carjhrj 

(7) 


Bab = (a aCefbr, 

(8) 


fX T'qri 

(9) 


Qa T-qa 

(10) 


Uab = Tab 

(11) 

Here eabc = ^rjabc where eajS'yS is the Minkowski spacetime volume element. Then 

eqn. (jHI) 

yields two constraint equations 




o'" Eab = 47ra(lcla(2/i -h 17%) 

- driQa) 

(12) 


d^Bab = 4:7iae''^ adeQf 

(13) 

and two equations of motion 




drjEab 

\ea^^daB,, - leb^^daB,a 


= dvra 

d(aqb) - ISabdcQ^ - d,{Uab - ^SabU^) 

+ 4na'{Uab - l6abU%) 

(14) 

d,Bab + \ea^^d,Eab + \eb"^d,E^ = 2na 

(15) 


Here a' = da/drj and 5ab^ the Kronecker delta, is the spatial metric of Minkowski spacetime. 

We now want to decompose the spatial tensors into tensors on the two-sphere. We intro¬ 
duce the usual spherical polar coordinates (r, 6*, 0) with the usual relation to the Cartesian 
coordinates {x,y,z). We use capital latin letters to denote two-sphere components. From 
the electric part of the Weyl tensor Eat we obtain a scalar as well as a vector and a 
symmetric, trace-free tensor given by 


— Eat ( 16 ) 

Eab = Eab — \HabEc^ ( 17 ) 

Here Hab is the metric on the unit two-sphere, and all two-sphere indicies are raised and 
lowered with this metric. Similarly, the decomposition of the magnetic part of the Weyl 
tensor yields Brr and 


hd — Bat 
Bab = Bab — ^HabBc^ 
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( 18 ) 

( 19 ) 






The decomposition of the spatial vector qa yields a two-sphere scalar and vector qa, while 
the decomposition of the spatial tensor Uab yields two-sphere scalars Urr and N = 
vector Va = Uav and a symmetric trace-free tensor 

Wab = Uab - \HabUc^ ( 20 ) 

Then the constraint equations (eqns. flT^ and flT^ ) become 

drErr + 3r~^E^r + r~^D^XA = dvra (^^dr{2qi + N) — drjQr'j (21) 

drBrr + 3r~^Brr + r~^D^YA = 4:Trar~^e'^^ DaQb ( 22 ) 

dr^A + 2r ^Xa — ^DaEj-t + r ^Eab 

= 47ra(iT>A(2/i + N) - d^^qA) (23) 

drYA + 2r~^YA - \DABrr + r~'^D^BAB = ^'KaeA^iDsqr - drqB) (24) 

Here Da is the derivative operator and cab is the volume element of the unit two-sphere. 
The evolution equations (eqns. (ITT)) and flT^ l become 

drjBrr + DaXb = 4:Tiar~'^e^^ DaVb (25) 

d^Err - T-^^^DaYb 

= Aira (drQr — dr/Urr + “ /W)) -|- Awa' (Urr — |(2iV -|- /i)^ (26) 

dr,YA + \r-h^^DcEBA + \eA^{3DcErr - 2drXc) 

= 27ra (eA^{\Dc{3U„ - N) - drVc) + r-h^^DBWcA) (27) 

d,XA - \r-h^^DcBDA - \eA^{3DcBrr - 2drYc) 

= 2-Ka{DAqr + drqA) - Ana{r~^qA + d^VA) + dvra'VA (28) 

d-qBAB + {DcXb + t' ^Ecb — drEcB) 

+ \eB^{DcXA + r-^EcA - 5.^ca) + ^HABe^^DcXo 
= 27iaeA^{D cVb + r ^Wcb — drWcB) 
+27iaeB^{DcVA + r ^Wca — drWcA) + 2'KaHAB^'^^D cVd (29) 

dqEAB — \^A^ {DcYb — drBcB + f ^Bcb) 

-\eB^{DcYA - drBcA + t-^Bca) - iHABe^^DcYn 
= dvra (^D^AqB) - \HABDcq^) - O^Wab) + dTra'lTAB (30) 
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CALCULATION OF MEMORY 


We now consider the behavior of the fields at large distances from the sonrce. Unlike 
the asymptotically flat case, we cannot make use of the formal dehnition of null inhnity: 
de Sitter conformal inhnity is spacelike, and all gravitational radiation is negligible there. 
Instead we dehne the optical scalar u = rj — r and consider the case of large r and moderate 
values of u. Note that in this case “large r” means large compared to the wavelength of the 
gravitational waves emitted by the source, but not large compared to the Hubble length. 
That is, we treat the case where ra' is of order 1. In the case of Minkowski spacetime, 
it is shown in [n| that stress-energy gets to large r and moderate u by traveling in null 
directions: that is, the dominant component of the stress-energy takes the form 


= AdaudfjU (31) 

We will assume that eqn. fl3T|) also holds in our case. From conservation of stress-energy, 
it then follows that A takes the form A = La~‘^r~^ where L is a function of u and the 
two-sphere coordinates. In physical terms, the quantity L is the power radiated per unit 
solid angle. That is we have 

= N = Urr = —Qr = La~^r~^ + ... (32) 


with all other components of the stress-energy falling off more rapidly. Here ... means 
“terms higher order in We assume that, as in the asymptotically flat case, the electric 
and magnetic parts of the Weyl tensor behave as follows: 


Eab = GABr + ... 

(33) 

Bab = i>AB'^ + ... 

(34) 

Xa = XAr~^ + ... 

(35) 

Ya = yAr~^ +... 

(36) 

Err = Pr~^ + . . . 

(37) 

Brr = Ql" ^ . 

(38) 


Here the coefficient tensor helds are functions of u and the two-sphere coordinates. Fur¬ 
thermore, we assume that, as in the asymptotically flat case, in the limit as |m| —>■ oo the 
only one of these coefficient tensor helds that does not vanish is P. Note that because of 
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the relation between Cartesian and spherical coordinates Eab behaving like r corresponds 
to Cartesian components of the electric part of the Weyl tensor behaving like r~^. 

Now keeping only the dominant terms in eqns. fl2HI24p and using the fact that ra' is of 
order unity, we obtain 


—SiiLa ^(a + ra') 

(39) 

—Q + D^ua = 0 

(40) 

-Xa + D^cab = 0 

(41) 

—ijA + D^Bab = 0 

(42) 


Here an overdot means derivative with respect to u. Similarly, keeping only the dominant 
terms in eqns. (I25tl3np yields 


Q + e^^DAXB = 0 

(43) 

P — e'^^DAUB = 87rLa“^(a + ra') 

(44) 

Va + \e^^DceDA + l^-A^xc = 0 

(45) 

Xa — DcBda — ^ga^Vc = 0 

(46) 

Bab + ^A^ecB = 0 

(47) 

Gab — ^A^bcB = 0 

(48) 


From here, the analysis proceeds essentially as in [h|. By convention, the scale factor a 
is unity at the present time. Therefore at the position of the detector, r] is the same as the 


usual time. Eat (despite the factor of a~^ in eqn. ((71)) is equal to the physical electric part of 
the Weyl tensor and thus is directly related to tidal force, and a' is equal to Hq, the Hubble 
constant. Note that eqn. fHHD is redundant, since it is equivalent to eqn. fH7p . Since cab 
and bAB vanish as u —)■ —oo, it follows from eqn. (I47ll that Bab = —^a^^cb- This can be 
used to eliminate hAB from eqns. fH2l) and fH6l) which then become 


Va + ^^^DcCda = 0 (49) 

xa - \D^ecA - \(^A^yc = 0 (50) 

Combining eqn. with eqn. fj45l) then yields 

Va + (^a^xb = 0 (51) 
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However, since xa and ha vanish as m ^ — cx), it then follows from eqn. (I5T|) that 


Va = -eA^XB 

(52) 

Thus, we can eliminate i/a from eqns. fl4UI) and 0441) which then become 


Q + e^^DAXB = 0 

(53) 

P — D^xa = 87rL(l + tHq) 

(54) 

But these equations are then redundant, since they are equivalent to eqns. 

flT3|) and fl3^ re- 

spectively. Thus the only independent quantities are xa, P, Q and L. 

satisfy the following equations 

These quantities 

D^cab = Xa 

(55) 

e^^DBCcA = (^A^Xc 

(56) 

Dax^ = P — 87rL(l + rHo) 

(57) 

e^^DAXB = -Q 

(58) 

Now let’s consider how to use eqns. 0551I58P to hnd the memory. Recall that cab is (up 

to a factor involving the distance and the initial separation) the second time derivative of 

the separation of the masses. Thus we want to integrate cab twice with respect to u. Dehne 

the velocity tensor vab, memory tensor uiab and a tensor 2 ;^ by 


pu 

Vab = / eABdu 

J —00 

(59) 

roo 

mAB = / VABdu 

J —CO 

(60) 

poo 

Za= XAdu 

J —CO 

(61) 


Now consider two masses in free fall whose initial separation is d in the B direction. Then 
after the wave has passed they will have an additional separation. Call the component of 
that additional separation in the A direction Ad. Then it follows from the geodesic deviation 
equation that 


Ad = — m^B 
r 

(62) 

To hnd rriAB we hrst integrate eqns. fl5^ and fl56|) to obtain 


D^vab = Xa 

(63) 

e^^DBVcA = ^A^Xc 

(64) 


Then integrating again from — cxo to cxo we obtain 


D^rriAB = (65) 

e^^DBmcA = ( 66 ) 

Now integrating eqns. (IFT)) and fl58|) from —oo to oo yields 

Daz^ = AP- 87rF(l + rHo) (67) 

e^^DAZB = 0 ( 68 ) 


where the quantities AP and F are dehned by AP = P(oo) — P(—oo) and F = J^^Ldu. 
In physical terms, P is the amount of energy radiated per unit solid angle. In deriving eqn. 
fl68|) we have used the fact that Q vanishes in the limit as |m| —)■ oo. Since is curl-free, 
there must be a scalar $ such that 2;^ = Da^- Then using eqns. fl67)l and fl6^ we hnd 


AP — 87rP(l -1- tHq) 

(69) 

D^ruAB = Da^ 

(70) 


DISCUSSION 

We now consider the physical implications of these results, and in particular of eqns. fl69l - 
[TOD. As in the asymptotically flat case, there are two kinds of gravitational wave memory: 
an ordinary memory due to sources that do not get out to inhnity and a null memory due 
to sources that do get out to inhnity. The ordinary memory is sourced by AP, that is the 
change in the radial component of the electric part of the Weyl tensor. The null memory 
is sourced by P, the energy per unit solid angle radiated to inhnity. However, in contrast 
to the asymptotically hat case, there is a factor of 1 -|- tHq multiplying P. Note that in 
cosmology, the wavelength of light from distant sources is redshifted by a factor of 1 -1- 2 ; and 
that to hrst order in z we have 1 + z = 1 + tHq. Thus, expressed in terms of P and r it 
seems that the null memory is enhanced by a factor of 1 However, r is not a directly 
observed property of a distant object: instead we observe the luminosity of the object and 
infer a luminosity distance di related to r by = r(l -1- z)r. Since the observed memory 
is given by eqn. fl62ll which has a factor of r~^, it follows that when expressed in terms of 
P and di, the null memory is enhanced by a factor of {1 + z)^. However, P is the energy 
radiated per unit solid angle as measured by the observer, who is at cosmological distance 
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from the source. Instead, one might want to calculate the local F (which we will call F\oc) 
as measured by an obsever who is sufficiently far from the source to be in its wave zone but 
still at a distance small compared to the Hubble radius. That is, F\oc is the F that source 
would have in Minkowski spacetime. Because the energy of radiation is diminished by a 
factor of 1/(1 + z), it follows that F = Fioc/(l + z). Therefore, when expressed in terms of 
Floe and the luminosity distance, the memory is enhanced by a factor of 1 + z. 

Finally, we consider the possible generalization of our memory calculation to the case 
where the background cosmology contains dust in addition to the cosmological constant. 
What makes eqn. ([3]) so simple in a de Sitter background is that the terms on the right 
hand side involve derivatives of the stress-energy, and that these derivatives vanish in de 
Sitter spacetime. However, as shown in 

of motion for the Weyl tensor involve the shear of the fluid four-velocity, and the equation 
for the Weyl tensor must be supplemented with equations for the shear. Thus one obtains 
a set of coupled equations involving the Weyl tensor and the shear. Nonetheless, it is 


argued by Thorne 


18| that for gravitational waves from realistic sources propogating in 


a cosmological spacetime, one can apply the geometric optics approximation in which the 
wavelength of the gravitational waves is small compared to all other length scales in the 
problem. The implication of this argument for the calculation of gravitational wave memory 
in the cosmological setting is that one expects all the extra terms in the equations due to 
shear to be negligible. This question is currently under study. 
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